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Abstract 



Using the QCD operator product expansion, we derive the real part of the 
transverse and longitudinal vector-vector correlation function with the quan- 
tum numbers of the p and u mesons to leading order in density and three 
momentum (q 2 ) for lu 2 — ► — oo. The operator product expansion provides, 
through the Borel transformed energy dispersion relation, a model indepen- 
dent constraint for the momentum dependence of the vector meson spectral 
density in nuclear matter. Existing model calculations for the dispersion ef- 
fect of the p, where the vector-meson nucleon scattering amplitude is obtained 
by resonance saturation in the s-channel, in general violate this constraint. 
We trace this to an inconsistent choice for the form factor of the ANp ver- 
tex. With a consistent choice, where both the form factor and the coupling 
constant are obtained from the Bonn potential, the contribution of the A is 
substantially reduced and we find good agreement with the constraint equa- 
tion. We briefly comment on the implications of our result for attempts to 
interpret the enhancement of low-mass dileptons in heavy-ion collisions. 
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I. INTRODUCTION 



It is expected that the properties of vector mesons at finite density and/or temperature 
should be reflected in the dilepton spectrum observed in nuclear reactions. It has also been 
suggested that modifications of the vector-meson properties are related to the restoration 
of chiral symmetry in dense matter 0,0- Indeed the enhancement of low- mass dileptons 
found in heavy-ion collisions H can be interpreted in terms of a non-trivial change of the 
vector meson spectral density in hot/dense matter. However, the dilepton spectra observed 
in heavy-ion collisions involve dileptons from different stages of the collision. Consequently 
collision models have to be invoked to relate the spectra to the vector-meson spectral density 
at a given temperature and density. At present, the data can be qualitatively interpreted in 
terms of two scenarios; on the one hand, a scenario where the vector meson masses decrease in 
matter and on the other hand, a scenario based on hadronic models where the p-meson 
spectrum exhibits a non-trivial dispersion effect (three momentum dependence) [|H],|n[] and 



an enhanced width 12,111 



In spite of attempts to unify the theoretical descriptions of the vector meson spectral 
densities in matter WM, the characteristics of the two approaches still differ and the rela- 



tion between them remains nebulous. In order to reduce the discrepancies and eventually 
approach a unified description, it is important to reduce the uncertainties of the theoretical 
models by constraining the parameters with experimental data [|l^,|l^] and results derived 
directly from QCD. 

In this work, we derive in detail QCD constraints on the three-momentum dependence 
of the vector meson spectral density in nuclear matter. These constraints have been derived 
previously in the context of QCD sum rules [|TT|,|T8! and shown to imply a small three- 



momentum dependence of the vector meson mass at finite density. Here, we give a detailed 
derivation of the constraints and test to what extent it is satisfied by existing hadronic 
model calculations. In the first stage we derive the Operator Product Expansion (OPE) 
and constructs constraints on the spectral density using the Borel transformed dispersion 
relation. In the second stage, model calculations are confronted with these constraints. As 
we will see, the relevant condensates that appear in the constraints for the dispersion effect 
are well known up to nuclear matter density and the corresponding QCD sum rule provides 
a solid reference point for model calculations ]P7|-|ID|]. 

An alternative approach is to model the phenomenological side by a simple functional 
form with a few parameters and obtain medium modifications of these parameters and thus 
of the spectral density. This procedure is in general not unique unless the form of the 
spectral density can be constrained by independent experimental or theoretical input. For 
the p meson, for instance, there is a trade-off between a reduction of the mass and an 
increase of the width. Thus, there is a band of possible combinations of mass and width, 
where the quality of the fit to the constraint is almost unchanged. Furthermore, there are 
uncertainties related to the medium expectation value of the four-quark operator, for which 
the ground state saturation hypothesis was used. Consequently, the resulting constraint || 
can be satisfied either by a large decrease in mass and a small increase in width or by a 
small decrease in mass and a large increase in width or something in between [|T2| , pO| . 



The paper is organized as follows. In section II, we derive the OPE relevant for the 
dispersion effect of the vector meson spectral density. In section III, we state the constraint, 
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which should be satisfied by model calculations. Finally in section IV, we discuss to what 
extent the existing hadronic model calculations satisfy this constraint. 



II. OPERATOR PRODUCT EXPANSION 



Consider the vector-current correlation function in nuclear matter 
n^(u;,q) = i J d'xe i(3 - x (G\T[J,(x)J u (0)]\G) , 



(1) 



where J M = qj^q, \G) is the nuclear ground state at rest and q^ = (uj,q). In what follows, 
when we give results for explicit vector mesons, we will use the currents = \{u^^w^d^^d) 
for the p and uj mesons, respectively. 

In general, because the vector current is conserved, the polarization tensor in eq. ([I]) has 



only two independent functions |^T 

n^(^, q ) = n T g^ + n^p 



[IV 



where we assume the nuclear medium to be at rest, such that P^ = 



Sij - q 4 qj/q 2 and V L = (q^q v /q 2 - - P 



polarization functions are given by 



fiu/ 



(2) 
: 0, 

The longitudinal and transverse 



P T 



hr = — n, 



00 



nr = -i(in;j + iiiflo) 



(3) 

q* I q* q* 

In the limit of q — > 0, the longitudinal and transverse polarizations are degenerate, i.e., 

n L = n T . 

In this work, we address only the leading three momentum (q) dependence of the po- 
larization functions. Consequently, we expand the correlation function in a Taylor series 
around q = and examine the dispersion relation at fixed q, 



Ren LiT (a; 2 , q 2 ) = Re (ll (uj 2 , 0) + IT], T (u 2 , 0) q 2 + 



du 



'p°lA u > q ) , pIA u >°) 



+ 



LO 



'v? — LO 2 ] 



q 



(4) 



where p(u, q) = (l/7r)ImIT R ('u, q), and R denotes the retarded correlation function. The 
OPE for 11° appropriately modified for the nuclear medium together with the corresponding 
energy dispersion relation provides constraints on the mass shift and spectral changes at 
q = H. In this work, we will explore the three momentum dependence by studying the 



OPE in medium for IT^ T . 



In general, the OPE f22j]23| for the polarization function at large Q 2 
the form, 



-uj 2 + q 2 is of 



1 



-Com|Q 2 |+£— A- 



t=2 k=l 
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2 2k 

X Qfj, 3 QfJ, 2 k Q4k+r-2 An-M2ft 

2 2fc 

+ E E " %1vl %2 k nik+T -2 C ll^2 k ■ ( 5 ) 

T=2 fc=l V 

Here we have extracted the trivial dependence so that A d ' T , C d,T represent the residual 
Wilson coefficient times matrix element of an operator of dimension d and twist r = d — s, 
where s = 2k is the number of spin indices of the operator. The first set of terms is due to 
the OPE of scalar operators, while the second set originates from operators with spin, which 
have been written as a double sum in twist and spin. The coefficients A, C represent the 
two linearly independent sums of operators, which reflect the two polarization directions. 

The longitudinal and transverse polarization functions can be obtained by using eq. ([5]) 
in eq. The q dependence emerging from the first line of eq. (|5[), namely the contributions 
from the scalar operators, is hidden in Q 2 . This we call the "trivial" q dependence, where the 
momentum q appears in the Lorenz invariant combination u 2 — q 2 . Here, we are interested 
only in the "non-trivial" q dependence of n^ T . Consequently the scalar operators do not 
contribute to T and only operators with spin contribute to the non-trivial q dependence. 
This is done by expressing the polarization function in terms of Q 2 and q 2 , i.e., n(Q 2 ,q 2 ) 
and extracting the term linear in q 2 . 



A. Linear density approximation 

To leading order in density 

(G\A\G) = (o\A\0) + £±<p\A\p) , (6) 

where \G) denotes the ground state of nuclear matter. The first term on the right hand side 
is the vacuum expectation value, which vanishes for operators with spin. In the second term 
\p) denotes a nucleon state with the normalization (p\\p') = 2p 6 3 (p — p'), p n the nuclear 
density and m the nucleon mass. As in vacuum, we include operators up to dimension 6 
in the OPE. This implies that in eq. (||) there are contributions from (r, s) = (2,2), (2,4) 
and (4,2). The nucleon matrix elements of the r = 2 operators are very well determined. 
The t = 4 matrix element appearing in the p, to sum rule are similar to those appearing in 



deep inelastic scattering (DIS) of electrons off nucleons [24]]. These matrix elements have 
been estimated f2*5|]2l| with an uncertainty of about ± 30% using DIS data from CERN and 
SLAC. 

The operators with spin indices are symmetric and traceless so that when the matrix 
elements are taken, the general form can be written as follows PTLES], 



M2 9^2j-it^2jP^2j+i ' ' ' P^2k ~\- permutations] J A2k , (7) 



\j=o 



where there are ( 2 fc— 2i)!2Jj! ^ erms f° r a given value of k and j and aj = (— 1) J ^ 2 k)\ ( more 
details are given in Appendix A). For s = 2 and s = 4 one finds 
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A 



liua/3 



P 



PfiPuPaP/3 ~ -^(PtiPu9af3 + PiiPa9u(3 + PuPpQocu + PvPa.9^ + PvPp9^a + PaP^9^) 



P t 

+ -^{9/iu9a/3 + 9iia9v(i + 9ii(i9v 



A, 



(8) 



Substituting eq. (0) into eq. (|^), we obtain the contributions of the spin dependent operators 
to linear order in density 



n =*l 

2m 



°UEE 



4 



2fc+r,r fc 
2k 



r=2 k~- 



i q 



r-2 



j=0 



Q 2 / a; 2fc - 2 J (2fc)(2fc - l)(2k - 2j - 2)!j! 



Cir T ' T J\fm 2 Y 1 (2fc-j) 



j,2fe+T,T fe / 2U 

+ ^ S S ~ J Q^ =r S lo^J * 2fc ^' (2* - 2j) !j! 



— e 



EE 



/( 2fc+T,T fc 

^2/c 



§ (<5 2 J s 2fc - 2j (2fe)(2A; - l)(2fe - 2j)!(j - 1)! 



2\ J 



1 



4-(2A;-j)-(2A;-j) 



r=2 fe= 



T-2 



(9) 



where x = = g^ v -q^q v /q 2 and = -p^Puq 2 /{p-q) 2 + {p^qv+Pvq^)/p-q-g^- The 

latter two are the usual polarization tensors used in DIS |23j. The sum over k corresponds 
to a sum over the number of spin indices while the sum over j is a sum over the target 
mass corrections. It should be noted that in the Bjorken limit x —>■ finite , Q 2 — > oo, the 
leading terms are twist-2 with any number of spins, while the higher twist terms and the 
target mass corrections are suppressed by powers of 1/Q 2 - 

The nontrivial contribution to H\ T can be obtained as follows. First, we use eq. (^) in 
eq. fl3|). Since we have assumed the ground state to be at rest, the nucleons are also at rest 



P 



(m, 0) so that eoo 



i2, doo 



2. We then extract the term 



proportional to q as discussed above. For the s = 2 terms we find 



ni = ^V(-i) r/2 ^— (4C 2 2 , . 



m 



l2+T,T\ 



T 2m V u T+i 2 2 1 



(10) 



while for s — 4, 



n! 



^(_ l) r/2+l 

2mV 



^(4<--20Q 4+ -) 



Hi = i2LV(_iWW 



m 



T+6 



;i84 +r ' T 



20C 4 4+r - T ) . 



(11) 



In general the contribution of operators with s spin indices is a sum of terms of order 
1/uj t+s+2 . 
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B. Form of A, C for r = 2 in the linear density approximation 

The coefficients A and C can be determined by comparing eq. (|S]), in the Bjorken limit, 
with the expectation value of the polarization tensor (Ef) in a one-nucleon state [p5j,p9fl , 



x' 



/ / d A xe^(p\T[J^x)J v (0)}\p) 

1 

Q 7 



(12) 



where in the MS scheme, the Wilson coefficients are given in Table dH). For a proton state 
normalized as (p\p) = (27r) 3 2po<5 3 (0), we have 



AUfi) = 2 fdx a^fafo/x) + {-l) n q{x^)\ 
Jo 

A G U) = 2 I dx x n ~ x G{x,n), 
Jo 



(13) 



where q(x, fi) (G(x,/i)) is the quark (gluon) distribution function inside the proton at the 
scale fi 2 . Comparing eq. (^) multiplied by |^ with the leading term in eq. © in the Bjorken 
limit, we find, 



4 2fc+2,2 
A 2k 



( - / 2,2fc /1 2fc ' L/ 2,2fc /1 2fc 5 



^2^+2,2 _ r q A q , n G aG 
°2fc — ^L,2k A 2k ^L,2k A 2k • 



(14) 



C. Form of A, B for r = 4 in the linear density approximation 



The contributions due to r = 4 operators have been computed in |24j and those due to 
the mass dependent r = 4 operator in 



z / d 4 xe^{p\T[J^x)J u (0)]\p) 



X 2 Q 2 



13 



d^H 1 + -H 2 + — H 9 - —H rn ) + e 



16 



4 



-H 2 



-H 9 - -H r 

1 2 



(15) 



where, for a current of J® q = qQqj^q, the H 's are the spin independent parts (such as the 
A in eq. @) of the following operators, 

H lu = {9 2 {qi^i^Q q q){qivi^—Q q q)) , 



fj,!/ 



" q=u,d,s ~~ 



H^^(ig(q{D^G va }j a l5 Qlq) 
H™ = (qD,D u m q Q 2 q q) . 



(16) 



Comparing eq. ((L5|) multiplied by I 2 - with the s = 2, r = 4 part of eq. @, we find, 
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H 1 + -H 2 + —H" - -H r 
8 16 4 

L H 2 - -H° - -H m . 
4 8 2 



(17) 



The matrix elements of these operators are not all known. However, the values of two 
different combinations in the F 2 and Fl directions have been obtained by one of us p5[ by 
analyzing the DIS data at CERN and SLAC. In order to proceed, we make the following 



assumption 25 



(18) 



where is some gamma matrix and A v is an isospin singlet operator and we take j3 = 0.48. 
Now, it is possible to uniquely determine the matrix elements appearing in the p, to meson 
sum rules. To see this, let us rewrite the scalar part of the matrix elements in eq. ([IB]) , 



H 1 = Q 2 u Kl + Q 2 Kl - (Q u - Q d ) 2 KlJ2 
H 2 = Q 2 u K 2 u + Q 2 d K 2 d , 
H° = Q 2 u Kl + Q 2 d Kl , 



(19) 



where we consider only the u,d sector and neglect the quark masses m u ,md- The K's are 
the matrix elements obtained from eq. ([15]) . 



Kl = —(p\ur + A\u\p), 



1,2 



rn 

U 9 



(p\u{D + ,*G + ,}^ l5 u\p) 



K d = —Ap\2{uT\u){dT\d)\p) . 



m 



(20) 



Here, r£ = 7/l75 f , r 2 = 7jU f , T + = (V + T 3 )/V2 and = uT^u + dT^d is a flavor singlet 
operator. 

Using the DIS data for Fl, F 2 on proton and neutron targets, one can now determine 
the three combinations 



Zi = K l + -K 2 + —K 9 ~ 

i u 1 g u ' \Q 



-0.06 GeV 2 



Z 2 = K 2 U - -K^ 0.56 GeV 2 



Z 3 = KlJ2 ~ -0.04 GeV 2 . 
The operator product expansion for y is then given by 



(21) 



^ • Pr. 

CO 8 



](1 + P)[-2Z 1 + ~Z 2 



1 ± 1 



(22) 



In the last term in Hl,t, the + sign refers to the p and the — sign to the 10 meson. 
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D. Final form of the OPE 



The final form of the OPE for II 1 can now be obtained by using eqs. (|T(], |TT], [14]) and 



n L,r(^)/Pn - 



L,T 



+ 



<jJ c 



For the longitudinal (L) and transverse (T) parts one finds fl~7|JT8 
1 



hi 



u 2 



Q 5 3 d 9 G G 3 G 

v^Cl,4 — -jCl^m A\ + (2^2,4 — 5C L4 )m A 4 



1 

+ -m 
2 



(l + /3)(-Z 1 + ^ 2 ) + (l±l)Z 3 



-CI 2 mA^ +d - C?,mA 



G 

L ) 2" b ^ L 2 1 



bl = {\CI A - \c%,)m*At d + (C& - SCgjm 3 ^ + ^(1 + (3)Z 2 , 
where again the ± refers to the p and u mesons, respectively. 



(23) 



(24) 



E. Infinite sum of twist-2 contribution 

Since the matrix elements of all the r = 2 operators are known, one can sum the higher 
order -^( n > 6) twist-2 correction. Here we present the result to lowest order in a s , i.e., 
including only the quark part of eq. (0). By substituting this into eq. (||) and evaluating 
the sums we find, 



m 



2oo 6 
2uo e 



p n m 



dxx- — — 
(1- 



m 2 x 2 \ 2 



[q(x) + q{x)} 



dx 



(1 + 



m 2 x 2 \(j 



(I" 



,,,2 



- 12 



(1 + 



(1- 



+ 3 



(1 + 



(1- 



m 2 x 2 \2 



x / dy[q(y) + q(y)} 



dxx 



(l + 22?£)6 (1 + 



m 2 x 2 \Q 



(I" 



m 2 x 2 \ /j 



+ 



(1 + 



(I" 



m 2 x 2 \2 



1 

x / dyy ~ x [q(y) + q(y)] , 



p n m r 1 (1 + 

— ^- / dxx 

uj 6 Jo (1 - 



m 2 x 2 \ 2 



[q(x) + q(x)] 



p n m 



dx 



16 



(1 + 



m 2 x 2 \5 



(1 



14 



(1 + 



(1 



m 2 x 2 \ 4 



+ 



(1 + 



(1 



m 2 a .2 ^ 



8 



p n m 



LL- 



16 



m 2 x 2 \fi 



x I dy[q(y) + q 



) 3 A 1 + 

— + 4- 



x / dyy [q(y) + q(y)} . 



(25) 

In the following section we use this partial summation to estimate the effect of higher 
dimensional (n > 6) operators on the constraint. 



III. CONSTRAINTS 

We now make a Borel transformation of uo 2 times eq. (^) and examine the dispersion 
relation, 

Pn {W - 2l^) = / ^ pV) ^ e ~" 2 ' M2 ~ Pnbscatt ■ (26) 

The reason for multiplying by u 2 is to get rid of any possible subtraction constants pro- 
portional to l/u> 2 . There could be further subtraction constants proportional to Pnbscatt/u 4 , 
which we include. In a given model calculation, the corresponding value of b scatt should be 
computed. The contribution from particle-hole intermediate states is b scatt = 1 /4m for the 
longitudinal p, u meson and zero for the transverse parts. (This is a non-trivial contribu- 
tion in the sense discussed above.) After the Borel transformation, the contribution from 
dimension=n operators to the OPE is divided by (n — 2)!. This improves the convergence 
of the expansion, and enhances the contribution of the low-mass region. 

In Fig. 1, we show the OPE side of eq. ( |2B| ) as a function of M 2 for the p meson. The 
solid lines shows the full calculation including dimension 6 operators, while the dot-dashed 
line shows the OPE without the twist-4 contributions and the long dashed line without 
dimension 6 contributions. At M^ in ~ 0.8 GeV 2 (M^ in ~ 2.5 GeV 2 ), the contributions of 
dimension 6 operators is approximately 40% of those from dimension 4 operators for trans- 
verse (longitudinal) polarization. Therefore, the OPE should be reliable only for values of 
the Borel mass above these, where one expects the contributions from higher dimensional 
operators to be suppressed. We also show (short-dashed line) the OPE including the con- 
tributions from the infinite sum of twist-2 operators, eq. (^). The corresponding Borel 
transform is given in Appendix B. The fact that the contribution of twist 2 operators of 
dimension larger than 6 is very small for M 2 > M^ in , indicates that the truncated OPE is 
indeed reliable above the minimum Borel mass. 

Thus, the constraints on the spectral density are given by eq. (|26|), for M 2 > M^ in . 
The Borel transformation also introduces an exponential weighting factor of the spectral 
density exp(— s/M 2 ), which for small values of Borel mass, suppresses the contributions 
at large energy. This has the advantage that in practical calculation, one can concentrate 
on the change of the spectral density near and below the vector meson mass region, while 
the contribution from higher energies can be effectively taken into account by a continuum. 
This approximation will be valid for Borel masses below some maximum value M max . We 
estimate M max by requiring that the continuum contribution should be less than that of the 
low energy region. 
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IV. COMPARISON WITH HADRONIC MODEL CALCULATIONS 



A possible non-trivial momentum dependence of the vector-meson self energy in nuclear 
matter has recently been extensively discussed. Such a momentum dependence arises natu- 
rally through the energy dependence of the p-nucleon scattering amplitude. In this context 
p-nucleon interactions in p-wave channels play an important role [|IO| , |rT| . |3~Tl . In a complemen- 
tary approach, the momentum dependence due to in-medium modifications of the p-meson 
pion cloud has been explored |3(J. We note that in none of these models is the momentum 
dependence directly constrained by data. Therefore it is useful to confront them with the 



constraints from QCD sum rules. Here we explore models of the type studied in refs. fllQjpl 
where the p-nucleon interaction is mediated by s-channel resonance exchanges. 

To make a systematic comparison with our constraints, we have to first relate the 
correlation operator eq. ([!]) and the polarization functions eq. (0) to the p meson self- 
energy T, LjT (ui,q), as computed in a hadronic model. To this end, we invoke vector-meson- 
dominance 

1 771^ ! 

n L>T (^, q) = -~^~ 2 ^— =— , r , (27) 

Q 9 P q-m 2 p + %Ym p - E Li t{uJ, q) 

where g p = 6.05 is the pirn coupling constant. 

The next step is to extract the lowest order nontrivial q 2 dependence to linear order 
in density. To this order, the self energy is of the following generic form for p-wave reso- 
nances lin, 

y _ 2 ( fRNp Y t w( M 2 S x E R (q)-m N 

ST " q VW) [F(q)] T'V-^q)-^ • (28) 

In this section we consider only transversely polarized vector mesons; in the longitudinal 
channel the expected momentum dependence is much weaker |T7| , pi| . For s-wave reso- 
nances 



3lj the leading q dependence in p8|) is replaced by a factor proportional to lu , 
^ 2 ( fRN P \ 2 rcV ms'S's E R (q)-m N 

Zt = uj \ [F(q)] — Pn — — -— - , (29) 



m p I 2 u 2 - (E R (q) - m N 



where the form factor is parametrized as 



A 2 

F(q) = — with A = 1.5 GeV , (30) 

A 2 + q 2 



and 



E R {q) = Vm| + q 2 - -r H . (31) 

Here mN,fn R are the masses of the nucleon and the resonances, and T R is the total width 
of the resonance and its momentum dependence will be modeled as in ref. [f[0[], 



t \ 21+1 

r fl (w, q) = r" ( (32) 
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In the phase-space factor, I = 0(1) for s- (p-)wave resonances and 



p = \J[(s — m 2 N — ml) 2 — 4m 2 v m 2 ]/4s 



PR = y[( m R ~ m N- m l) 2 ~ 4mX]/ 4m Ji , (33) 

where s = (uj + m^) 2 — q 2 . However, in the present context the effect of the momentum 
dependence of the width turns out to be weak. 

As can be seen from eq. fl2~8|), the contribution from p-wave resonances is proportional 
to q 2 . Hence its contribution to Imn 1 is trivially obtained by neglecting the remaining q 2 
dependence. In contrast, eq. ( p9|) is proportional to u 2 . The contribution of such terms to 
Imn 1 is most easily obtained by replacing u 2 by q 2 + q 2 and expanding to leading order in 
q 2 keeping q 2 fixed. Finally, we obtain 



-t 4 
I TYl 

imnV 2 )res = 



1 

X 



(u 2 — m 2 + inipTp) 2 



( fwfpV S s m R -m N - %T R 



; , P -wave \ m P J 2 uj 2 -(m R -m N - \ T R ) 2 
. , , 2 ^s 2 ( -M m R- m N - 1 V R 



j2 ( lm] 

-wave \ m p / 



2 rn \oo 2 - (m R -m N - |r R ) 2 

-^h( m N + %Tr)(™> r -m N - |r 



2m R . m R 



uj 2 - {m R -m N - ^Yr) 2 (to 2 - (m R - m N - ^T R ) 2 ) 2 




(34) 



where we have neglected the small contribution due to the momentum dependence of Tr. 

The contribution from the higher energy region is modeled using the following assump- 
tion. In the vacuum QCD sum rule approaches, the spectral density is approximated by a 
pole and a continuum. Since we make a hadronic calculation of the changes in the pole, 
we will add the simplest change in the continuum and allow a non-trivial change in the 
continuum threshold, 



cod (q 2 ~ (so + s'q 2 )) , (35) 

where s' is a parameter to be determined. The parameter cq, which characterizes the strength 
of the high-energy continuum, is determined by QCD duality. This excludes a simple q 2 
dependence, such as c — > c + c'q 2 . We thus find the following contribution to the spectral 
density, 

Imn 1 (o; 2 ) corii = -c s'5(u 2 - s ) (36) 



The phenomenological spectral density p 1 is then given by the sum of eq. ([34]) and eq. (j36|), 
i.e., p 1 = (l/n)(lmIll es + lmIll ont ). The parameter s' in eq. (|36|) is determined by minimizing 
the difference between the left and right hand side in eq. (|2~6"D . We note that the parameter 
s = 1.43 GeV 2 is fixed by the QCD sum rule § at q = 0. 
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A. A 'standard' set of model parameters 



We first explore the model employed by Peters et al. [31] . The corresponding parameters 
are given in Table 2. We choose this set, since it is representative for the sets used also by 



other groups, with similar values for the parameters (see e.g. refs. ||lT|Jl6| , |33| ) . In Fig. 2 we 
compare the OPE with the phenomenological side, i.e., the left and the right hand sides of 
eq. (|26|) , including both the resonance and continuum contributions in the phenomenological 
spectral density. The best fit was obtained in the Borel window 0.5 — 2.0 GeV 2 , by adjusting 
the parameter s' of the continuum. We find a huge discrepancy between the OPE and the 
phenomenological side. The major contributions to the phenomenological side are due to the 
nucleon as well as the A(1232), A(1520), A(1680), A(1720) and A(1905) resonances. For 
comparison we show the right hand side for only these six states (denoted as 'selected subset' 
in Fig. 2) as well as some of the individual contributions. Clearly, the very large effect of the 
A(1232) is mainly responsible for the disagreement. In the best fit, its contribution, which 
in itself is much larger than the OPE side, cannot be compensated by the continuum over a 
reasonable range of Borel masses. The fact that the continuum lies at a much higher energy 
than the A(1232), leads to a characteristic difference in the dependence on the Borel mass. 
Thus, in most of the Borel window, the best fit lies far outside the estimated boundary of 
uncertainty shown in Fig. 1. 

In Fig. 3, we show the corresponding spectral density p 1 . The contribution of a particular 
resonance is identified by extracting the corresponding term in the curly brackets in eq. 



Also here the dominance of the A(1232) is apparent. Moreover, in the constraint equation 



([26]), the contribution of the A(1232) is further enhanced compared to the other structures 



at higher energy due to the exponential factor. Thus, we conclude that the strong effective 



ANp coupling originally employed by Rapp et al. |TTJ] is incompatible with the constraint 
from QCD sum rules. 

B. A consistent model 

As noted above, the most important contributions to the sumrule are due to the selected 
set of resonances. Consequently, in the following analysis, we will include only these baryon 
states. Let us first recall how the coupling constants fnN P are fixed. The couplings for 
the A(1720) and the A(1905) are determined from their decays into a nucleon and a p 
meson [|K|. Hence, these couplings should be reliable for the process p + N — > resonance 
considered here. If the mass of the resonance lies below the p N threshold, however, the 
corresponding coupling constant can clearly not be determined in this way. Consequently, 
in many models the p N couplings to the nucleon and the A(1232), are obtained e.g. from 
the Bonn potential |34|. An important point, which apparently has been overlooked, is that 



in the Bonn potential the coupling constant Jrn p is multiplied by a form factor 

V A 2 + q 2 

where A p = 1.4GeV and n p — 1 for the NNp and n p = 2 for the ANp vertex. Hence, the 
effective value of the coupling constant at q = is reduced by 0.7 for the nucleon and by 



FM 2 ) = {^r^) • (37) 
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0.49 for the A(1232). Since the self energy is proportional to the square of the coupling 
constant and the corresponding form factor, the contribution of the nucleon (A(1232)) is 
reduced by a multiplicative factor of 0.49 (0.24). This reduction of the effective coupling 
constant has been left out in previous calculations. 

The resonances iV(1520) and iV(1680) are situated below the nominal threshold for the 
decay into a nucleon and a hypothetical zero-width p meson. Thus, they can decay into a 
nucleon and a p meson only because of the finite widths of the particles involved. Since in 
this case only the low- mass tail of the p meson is seen in the two-pion spectrum |35|| , a reliable 



identification of the p is very difficult. Consequently, the determination of the corresponding 
coupling constants is rather uncertain. In order to explore possible consequences of this 
uncertainty, we also show results for the case where the two coupling constants vanish. 

Let us now confront this set of coupling constants with the QCD sum rules. For the 
nucleon and the A(1232), we use the couplings of the Bonn potential, with the corresponding 
form factors, extrapolated into the time-like region 



if q 2 < mj. 

F p c (q 2 ) = { Kll ~ P f r / m 2 V > (38) 
1 + ( 9 A ™ p ) otherwise. 



p 



Thus, for q 2 < m 2 we use the covariant generalization of (|37|). However, close to the pole 



at q 2 = A 2 this form clearly makes no sense. For q 2 > m 2 we therefore employ the simple 



form shown in the second line of (|38|), with the exponent n» = 0, 1 or 2. Fortunately, due to 
the exponential weighting, the QCD sum rule (|26|) is insensitive to the precise form of F c 
in this range of energies; the real parts of II 1 obtained from the right-hand side of ( p6[ ) for 
different values of n, are almost indistinguishable. Thus, for the nucleon and the A(1232) 
we replace the coupling constant f RNp in eq. (|341) by f RNp (uJ 2 ) = F p c (u 2 )f RNp . 

For the A(1905) and the iV(1720) we use the coupling constants f R N P determined in 



ref. ||10|| . Since the corresponding form factors are normalized to unity for zero three- 
momentum, they do not affect the leading q 2 dependence. Furthermore, a possible frequency 
dependence can also be neglected, because, in contrast to the A(1232), these resonances 
contribute significantly only in the kinematical region where the coupling constants are de- 
termined, i.e., for almost on-shell p mesons. The coupling constants used are given in Table 
III. Note that here we use the NNp and A(1232)A r p coupling constants of ref. |j53| , which 
differ slightly from those used in ref. [|31]]. 

In Fig. 4 we show the effect of a consistent treatment of the ANp form factor on the 
spectral density. As expected, the contribution due to the A(1232) is substantially reduced. 
This reduction is also reflected in an improved agreement with the constraint equation, as 
shown in Fig. 5. There, the distance between the two solid lines, which correspond to the 
OPE with and without the dimension 6 operators, represent the uncertainty of the left-hand 
side of (^6]). The long-dashed line is the right-hand side without the form factor for the delta 
and the nucleon, while the dashed line shows the result including form factors of the form 
eq. ([!8|) with rii = 0. (Since, as noted above, the dependence on rii is minimal, we do not 
show results for rii = 1,2.) Furthermore, the dotted line shows the result of including the 
form factors but with vanishing couplings to the iV(1520) and iV(1680) resonances. Thus, 
the range between the dashed and dotted lines indicate the uncertainty due to the poorly 
known iV(1520) and iV(1680) coupling constants. 
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Using the criteria discussed in Sect. Ill, we obtain a Borel window, where the constraint 
equation can be applied, of 0.8 GeV 2 < M 2 < 1.8 GeV 2 . Clearly, there is a good agreement 
between the OPE and the phenomenological calculation, when the ANp form factor is 
consistently taken into account. We note that most of the disagreement is removed already 
with a static (frequency independent) form factor, i.e., by renormalizing the effective coupling 
strength by 0.7 for the nucleon and by 0.49 for the A. The resulting phenomenological 
side is shown by the dash-dotted line in Fig. 5. However, the frequency dependence, as 
parameterized in (|38|), leads to a further improvement (see Fig. 5). 

Our work indicates that the ANp coupling used in many model calculations is unphys- 
ical. The Bonn potential offers no justification for the very strong effective couplings used. 
Furthermore, such couplings violate the constraints obtained from QCD sum rules. Thus, it 
is necessary to re-analyze those calculations of the low-mass dilepton spectrum in heavy-ion 
collisions, where the ANp coupling was treated inconsistently. 



V. DISCUSSION 



We have presented a model independent constraint for the three momentum dependence 
of the vector meson spectral function in a nuclear medium. We have shown that the mo- 
mentum dependence of the spectral density obtained in model calculations |Tl|,p!6|,|3T|,|33[] is 



too strong to satisfy the constraint derived from QCD sum rules. The main source of the 
discrepancy is the large contribution from the A(1232). In these models, the coupling of 
the rho meson to the nucleon and the A(1232) is not accounted for in a consistent manner. 
While the coupling constant is taken e.g. from the Bonn potential, the form factor is not; the 
normalization of the form factor differs considerably from that, which is used in the Bonn 
potential. This leads to an effective coupling strength which is approximately a factor two 
too large. When the original Bonn form factor is employed, the contribution of the A(1232) 
is substantially reduced and we find good agreement with the constraint equation. 

In view of this it would be important to reexamine the hadronic model calculations that 
are used to interpret the low-mass dilepton enhancement found in relativistic heavy-ion 
collisions, taking the form factors properly into account. The unphysically large coupling to 
the A(1232) gives rise to a spurious enhancement of the vector meson spectral density for 
low invariant masses. 
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APPENDIX A: COEFFICIENTS OF TRACE TERMS 



Here we specify the trace terms in detail. As discussed in the text, 
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' ' ' 9^2j-i^2jP^2j+i ' ' ' Vn2k ~\~ Permutations] j A<ik , (Al) 
where for a given value of k and j there are p8| (9 , io fc |i 9J - t different terms and a 



v 23 (2fc)! ' 

Now consider the two indices \xi. They can enter either as indices of the metric tensor 
g^ v or of the nucleon four momentum p^. Depending on how the two indices appear, we can 
distinguish four types of terms 

9 mi* x [{j - l ) factors (g^)] x [2(k - j) factors (p M )] 
Quirf/taX x [(j ~ 2 ) factors (g^)] x [2(k - j) factors (p M )] 
0mi«Pa« x [(J - 1) factors (g^)] x [(2A; - 2j - 1) factors (p„)] 
PwPw x [(j) factors (gr^)] x [(2fc - 2j - 2) factors (p M )] . (A2) 

There are (2fc_2j)!2^"-i(j-i)! numDer °f terms of the first type, where the remaining in- 
dices are ordered differently Of the second type there are (2k-2%z>-2(j-2y. ' °^ ^ ne third 

(2fc-2j-i) ; !2J 2 - !l (j-i)! anc ^ °^ ^ e f° ur t n there are ( 2 k-2j-2)\23 numDer °f terms. Using these 
results, the infinite series of twist-2 contributions can be summed. 



APPENDIX B: BOREL TRANSFORMATION OF THE INFINITE SUM 



The Borel transformations of eq. (25) can be obtained from the basic formula 



B.T. 



(oo 2 - A) 1 



1 



1 



(a - 1)! (M 2 )^ 1 



-A/M 2 



We find, 

B.Tyn ) 



dxx 



1 -1 _ -m 2 x 2 /M 2 \ 1 _2_ -m 2 a; 2 /M 2 

j + M 2 



m 2 x 2 



[q(x) + q(x)} 



+ dx 
J 



m 2 x 2 



m A x 4 , 



3 M 6 



'-^ _ e ~m 2 x 2 /M 2 ^ _|_ e -m 2 x 2 /M 2 

+ g(y)] 

+ 



2 2 

m x 



M 2 12 "M 4 



2 2 o 4 4" 



^e—/M 2 ^ _ 12 ^ + ^ J y dw - 1[ff(y) + ^ 



:bi) 



:b2) 



B.T.(^ 2 nj 

p n m 



+ / dx 
h 



2M 2 
6 

m 2 x 2 

,6^.6 



m 2 x 2 



+ e 



-m 2 x 2 /M 2 



+ 



e _ m 2 x 2 /M 2 _ 



-m 2 x 2 /M 2 



2rr2 J^2- 



\q[x) + 



(- 



6 



- 16- 



4 4 
m x 



+ 



32 m b x b 32 m 8 x* 



3 M 8 15 M 10 ' 2M 2 



M 2 M 6 
dy[g(y) + g(y)] 
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f 1 - m 2 x 2 /M 2 ( 32 m V 32 m 6 a; 6 32m¥\ 

dyy 1 {q(y) + q{y)] , 



el 



+ 



1 

da; 



m 2 x 2 



l 

_|_ g— m^ar 



/M 2 ^ 



1 2 

+ 



m 2 x 2 M 2 



o 



m?x 2 
200 mV 
3 M 6 



^ _ g -m 2 :c 2 /M 2 ^ _|_ e -m 2 x 2 /M 

64 m 8 a; 8 



[g(ar) + q(x)] 
6 



+ 36 



2 2 

m x 



+ 32 



m 6 a; 6 



+ [ 1 dxxe- 2 * 2 ^ f-24^ 



15 M 10 

176 m 4 a; 4 



M 2 \I> 
dy[q(y) + 



+ 



M 4 3 M 6 

'.X' 



m 6 x 6 64 m 8 x 8 
-32^- + 



M 8 15 M 10 
dyy- l [q(y) + q(y)} . 
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TABLES 



TABLE I. Wilson coefficients for r = 2. T(R) = |, C 2 (R) = §. 



Wilson coefficient 


n = 2 


n = 4 


C{ n = 1 + ^ 
Cl, n = ^C 2 (R)^ 

nG _a»T^D^4[ 4 4 i 1 n 2 +n+2 M , s^n 1 \i 
°2,n 47r i K K )f[(n+l) (n+2) 1 n y nfn+l)^)! 1 + ^j=l j)J 
(~<G _ o s 7V p\ 1 16 
^L,n ~ 47T 1 \ H )f (n+l)(n+2) 


1 + ^(0.44) 

a s 16 
4?r 9 

^l(_I) 
47I-V 2/ 

a s 2 

4?r 3 


1 + ^(6-07) 

a s 16 
4?r 15 

a s ( 133 \ 
4tt ^ 180 / 
a s 4 
4?r 15 



TABLE II. The parameters used in ref. [31]. 







^Total{GeV) 


l P 


fRNp 




N(940) 


1(1+) 

2^2 > 





1 


7.7 


4 


N(1520) 


1(3-) 
2V2 > 


120 





7 


8/3 


N(1650) 


2^2 > 


150 





0.9 


4 


N(1680) 


l/5 + \ 
2^2 > 


130 


1 


6.3 


6/5 


N(1720) 


1(3 + ) 
2V2 > 


150 


1 


7.8 


8/3 


A(1232) 


3/3 + \ 
2\2 > 


120 


1 


15.3 


16/9 


A(1620) 


3(1-) 
2V2 / 


150 





2.5 


8/3 


A(1700) 


3 / 3 — \ 
2^2 > 


300 





5.0 


16/9 


A(1905) 


3/5 + \ 
2\2 > 


350 


1 


12.2 


4/5 



TABLE III. Parameters used in a consistent model for the selected subset of resonances. 



I(J P ) r Total (GeV) f RNp 

N(940) \{\ + ) 8?T 

N(1520) ±(§~) 120 7 

N(1680) 130 6.3 

N(1720) 150 7.2 

A(1232) |(| + ) 120 16 

A(1905) §(| + ) 350 9.0 
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FIG. 1. The OPE in GeV -2 . The solid line shows the full calculation, while the short-dashed 
line includes the infinite sum of the twist-2 operators. The dashed line shows the OPE without 
dim-6 and the dash-dotted without twist-4 operators. All curves are at nuclear matter density. 
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FIG. 2. Real part of n 1 in GeV -2 . The solid lines denote the OPE and the best fit including all 
resonances in the model of ref. [31]. The contributions from the nucleon, A(1720), A(1232), A(1905) 
and the continuum are also shown. Furthermore, the contribution of the most important resonances 
(the selected subset) is shown for the parameters given in Table II. 
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FIG. 3. Imaginary part of II 1 = p 1 in GeV~ 2 as a function of the energy ui. 



0.5 

p 1 0.0 



„- Without form factor 




A A 




; \ A f\ 








/ v 

With form factor 1 









co [GeV] 

FIG. 4. Imaginary part of II 1 = p 1 in GeV -2 for the selected subset (Table II) with and without 
form factors (38) for the nucleon and the A(1232). 
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FIG. 5. Real part of II 1 in GeV -2 with the parameters given in Table III. The long-dashed line 
shows the phenomenological side without the form factor for the nucleon and the delta, while the 
dashed line shows the result including the form factors (38). The dotted line is obtained with the 
form factors but with zero coupling constants for the iV(1520) and iV(1680) resonances. Finally, 
the dash-dotted line shows the phenomenological side obtained when the frequency dependence of 
the form factors (38) is neglected. 
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